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Abstract 

^"^ The study of short-duration gamma-ray bursts (GRBs) experienced a complete revolution in recent years thanks to the discovery of 

^-^ the first afterglows and host galaxies starting in May 2005. These observations demonstrated that short GRBs are cosmological in 

^N^ origin, reside in both star forming and elliptical galaxies, are not associated with supernovae, and span a wide isotropic-equivalent 

^ energy range of ~ 10"^^ - 10^^ erg. However, a fundamental question remains unanswered: What are the progenitors of short 

gj GRBs? The most popular theoretical model invokes the coalescence of compact object binaries with neutron star and/or black 

^ hole constituents. However, additional possibilities exist, including magnetars formed through prompt channels (massive star core- 

^^^ collapse) and delayed channels (binary white dwarf mergers, white dwarf accretion-induced collapse), or accretion-induced collapse 

of neutron stars. In this review I summarize our current knowledge of the galactic and sub-galactic environments of short GRBs, 

i~n and use these observations to draw inferences about the progenitor population. The most crucial results are: (i) some short GRBs 

W explode in dead elliptical galaxies; (ii) the majority of short GRBs occur in star forming galaxies; (iii) the star forming hosts of 

p^ short GRBs are distinct from those of long GRBs, and instead appear to be drawn from the general field galaxy population; (iv) the 

cH physical off'sets of short GRBs relative to their host galaxy centers are significantly larger than for long GRBs; (v) there is tentative 

Q^ evidence for large off'sets from short GRBs with optical afterglows and no coincident hosts; (vi) the observed off'set distribution 

' is in good agreement with predictions for NS-NS binary mergers; and (vii) short GRBs trace under-luminous locations within 

^ their hosts, but appear to be more closely correlated with the rest-frame optical light (old stars) than the UV light (young massive 

^ stars). Taken together, these observations suggest that short GRB progenitors belong to an old stellar population with a wide age 

C3 distribution, and generally track stellar mass. These results are fully consistent with NS-NS binary mergers and rule out a dominant 

population of prompt magnetars. However, a partial contribution from delayed magnetar formation or accretion-induced collapse 

fvq is also consistent with the data. 
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1. Introduction 

. Gamma-ray bursts (GRBs) are short, intense and non- 
repeating flashes of y-ray radiation originating at cosmologi- 
jcal distances. While GRBs exhibit a broad diversity in their 
prompt y-ray emission (e.g., duration, spectral shape, peak en- 
ergy, brightness), they can still be divided into two basic cat- 
egories: shor t-duration and long-dur ation with a separation at 



about 2 sec ( Kouveliotou et al.L I1993D . The short GRBs have 



durations as short as ~ 10 msec, while the long events extend 
to hundreds of seconds. In addition, short GRBs tend to exhibit 
harder y-ray spectra than the long-duration events, and gener- 
ally have a lower fluences. 

The basic bimodality of GRB durations provided an early 
clue that the progenitors of the two classes are likely to 
be distinct. Within the broad range of possible scenarios, 
two popular models have emerged: The collapse of rapidly- 
rotating massive stars ("collapsars"; dMacFadyen and WoosleyL 
ri999.) ) in the case of long GRBs, and the coalescence of 
compact object binaries (with neutron star and/or black hole 



constituen ts - NS-NS/NS-B H; dEichler et al.lll989l: IPaczvnski' 
[1991; Nar avanetal.lll992h ) in the case of short GRBs. The key 
attractions of this mapping are the potential for a large energy 
release from both progenitor classes, and the expected typical 
timescale for each progenitor: A free-fall timescale of tff ^ 
30 s (M/10 Mq)-^^^(R/10^^ cm)^/2 f^^ collapsars, and a dynam- 
ical timescale of milliseconds for the compact merger remnants 
of neutron stars and black holes. However, other progenitor 
systems have also been proposed for short GRBs, for example 
magnet ars, thought to be the pow er source behind soft y-ray re- 
peaters (IThompson and Duncanlil995). accretio n-induced col- 
lapse (AIC) of neutron stars (iOin et al.L Il998h . and delayed 



magnetar formatio n throug h binary whi te dwarf merge rs or 
white dwarf AIC ( Levanetal. . 2006; M etzger et al.ll2008k . 

Until a decade ago, the distances, energy scale, geome- 
try, environments, and progenitors of GRBs, as well as the 
relation between the two burst classes, remained uncertain 
due to the lack of precise positions. The discovery of long- 
wavelength, long-lived "afterglows" from l ong GRBs in 1997 



provided the first glimpse at these properties (ICosta et al.L 119971: 
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Frail et al.L 119971: Ivan Paradijs et al.L 119971) . Indeed, the sub 



arcsecond positions enabled by long GRB afterglow detections 
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demonstrated a cosmological origin (IMetzger et alJ, Il 997h. an 
energy s cale of ~ 10^^ er^ (JFrail et all |2001; Berg er et al 



L2003ab; Blo om et all l20Q3k, significan t collimation (~ 10° 
jets; (Har risoneTall I1999I: IStanek et all [1999.)), and direct ev- 



idence for relativistic expansion (ITaylor et al. 



2004 ). Intense 



afterglow and host galaxy observations also linked long GRBs 
with the deaths of massive stars, mainly t hrough their exclu- 
sive location in star forming galaxies (e.g., (IBloom et all 1 19981: 



IPjorgov ski et all ll998':'Fr uchter et all 1 1999k ). their strong cor- 
relation with t he res t -frame ultravio l et (UV ) light of their hosts 
(IBloom et all I2OO2I: iFruchter et al.[ I2OO6I) , and their associa- 
tion with Type Ic co re-collapse supernovae (IHjorth et alll2003l: 



Stanek et al 



2003) 



The afterglows of short GRBs were discovered only in 2005. 
Prior to that point only a few afterglow searches were possible 
due to the relative faintness of the y-ray emissio n and hence a 
low e vent rate and large and delayed error circles (IHurley et al.[ 
12002"). In retrospect, these searches were woefully inadequate, 
reaching only about 21 mag at ^^ ^ 1 day in the optical and 
~ 0.5 mJy at dt ^ few days in the radio. The launch of NASA's 
Swift satellite in late 2004, provided the first chance for rapid 
and accurate positions for s hort GRBs, and indeed led to the 
detec tion of th e first X-rav (iGehrels et all 120051: iBloom et al. . 
2006h. op tical ('F ox et al.1 120051: iHiorth et all l2005k . and radio 
(lBergeretal.. ,2005l) afterglows in 2005. 

As in the case of long GRBs, the determination of accu- 
rate positions revolutionized the study of short GRBs. First, 
it led to an association with galaxies at cosmological distances 
(e.g.. teerger et all l2005i IFox et all l2005l iHiorth eTall l2005l : 



Bloom et all I2OO6I: Prochaska e t alll2006l: lBerger et all 120071) ) 
and henc e an energy scale of ~ 10"^ ^ - 10^^ erg (assuming 
isotropy) dBergeii 120071: iNakail l2007h . Second, it led to the 



and detailed properties of their host galaxies, and their locations 
within their hosts. The structure of this review is as follows. In 
^I summarize the discovery of short GRB afterglows, and the 
subsequent identifications of their host galaxies and redshifts; 
the detailed properties of the hosts (luminosities, metallicities, 
star formation rates, masses, stellar population ages) are dis- 
cussed in ^ and iH in ^ I discuss the sub-galactic environ- 
ments of short GRBs, utilizing mainly high-resolution Hubble 
Space Telescope observations; I discuss the possibility of large 
progenitor off'sets (due to kicks or globular cluster origin) in ^ 
and finally, in ^I use these results to place constraints on the 
progenitor population. 

2. The Discovery of Short GRB Afterglows, Host Galaxies, 
and Redshifts 

The discovery of short GRB afterglows starting in May 2005 
led to the first identifications of their host galaxies and hence 
to distance measurements. The first short GRB with an after- 
glow detection, GRB 050509B, was localized to a positional 
accuracy of about 5 ' ' with the Swift X- ray Telescope (XRT) 
( GehrelsetaD. l2005l : iBloom et all l2006h . No optical or radio 
emission was detected. The X-ray error circle appeared to co- 
incide with the outskirts of an ell i ptical galaxy at z = 0.22 6 
( Gehrelsetall 120051: IHiorth et all I2OO5I: iBloom et all l2006h . 
with a probability of chance coincidence of ~ 10"^. However, 
the error circle contained additional fainter galaxies possibly at 
higher redshift. 

Only two months later, GRB 050709 was the first short burst 
localized to sub-arcsecond precision through the detec t ion of 



X-ray (with Chandra) and optical emission ( Fox et all 120051: 



association of some events with elliptical galaxies po inting to 
an old progenitor populat ion (IBerger et all 12005: Gehrels et al. , 
l2005l:lBloom et alll2QQ6l) . Third, it demonstrated that the after- 
glow emission is similar t o that of long GRB s, albeit with a 



long 



Hjorth et all |2005). The resulting afterglow position coin 
cided with the outski r ts of an irregular sta r forming galaxy at 
z = 0.161 dFox et all Eo05l: ICovino et all Eob6). Despite the 
on-going star formation activity within the host galaxy, the burst 
was not accompanied by a supernova explosion, indi cating that 



generally lower luminosity (JBerger et all|2005|). F inally, it pro- the p r ogenitor was not lik ely to be a massive star dFox et al 



vided a rough estimate of the short GRB event rate (INakar et al. 
12006 *). 

Despite these fundamental results the progenitors of short 
GRBs remain unidentified at the present. The key observational 
test of the NS-NS/NS-BH merger model, the detection of coin- 
cident gravitational waves, is at least several years away. Simi- 
larly, theoretical predictions of early optical/UV emission from 
an accompanying "mini- supernova" dLi and P aczynskii Il998h . 
caused by the ejection of radioactive material from the merg- 
ing system, are highly uncertain, and even the most optimistic 
predictions lead to a faint and rapidly-fading signal that may 
be challenging to detect. Thus, the most promising avenue for 
progress at the present comes from statistical studies of the en- 
vironments of short GRBs, both on gal actic and sub-galactic 
scales (e.g., dProchaska et all I2OO6I: iBerger, 2009; Fong et al., 
2OIOI) ). These studies benefit from many of the same techniques 



20051: iHjorth et all l2005l) . However, due to the presence of 



active star formation, an association with a young progenitor 
system such as a magnetar could not be excluded. 
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that linked long GRBs with the death of massive stars, and al- 
low for comparison with theoretical predictions. 

In this review, I present the current state of our knowledge 
about the redshift distribution of short GRBs, the demographics 



Figure 1: Discovery images of the near-infrared afterglow of GRB 050724 and 
its elliptical host galaxy. The inset shows the Very Large Array radio posi- 
tion (ellipse) and the Chandra X-ray position (circle). This was the first short 
burst to unambiguo usly establish a link with an old stellar population. From 
IBerger etal] ( l2005h . 

It was only 15 days later that the discovery of X-ray, op- 
tical, and for the first time radio afterglow emission finally 
established a direct link between a short GRB and an old 



stellar population (IBerger et all l2005l: iBarthelmy et all |2005|; 
Gorosabel et"alll2QQ6h . The afterglow of GRB 050724 was lo- 
calized to an elliptical galaxy at z = 0.257 with no evidence for 
star formation activity (< 0.05 Mp V^^" ^ ) ^^^ with a stellar pop' 



ulation age of > 1 Gyr (IBerger et all l2005l: IProchaska et al. 



L2OO6). The absence of both star formation activity and an as- 
sociated su pernova demonstrat ed a direct link to an old stellar 
population (IBerger et alll2005h . 
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Figure 2: Discovery images of the near-infrared afterglow of GRB 070724 and 
its star forming host galaxy. Each row shows an afterglow image, a subsequent 
template image, and a residual image. The fading afterglow coincides with the 
disk of the host galaxy. From lBerger et al.1 ( 120091) . 

The combination of low redshifts (z ~ 0.2) and the ap- 
parent dominance of elliptical galaxies in the first few short 
GRB hosts led to initial speculati on of a particularly old pro- 
genit or population: r > 4 Gy r ( Nakaret all l2006h . r > 1 
Gyr (Zheng and Ramirez- Ruizi 12007 ), and r > several Gyr 
(jGal-Yam et al., 2008) . Indeed, a possible inconsistency with 
the expect ed merger time dela y distribution of NS-NS binaries 
was noted (JNakar et al.Ll2006h , although subsequent population 
synthesis models of NS-NS binary forma tion and mergers led to 
opposite claims (IBelczynski et al.L I2OO6D . Clearly, the sample 
of short GRBs with afterglow detections available when these 
various claims were published was very small (GRBs 050509B, 
050709, 050724, and 05 1221 A). 

Fortunately, the continued detection of short GRBs (mainly, 
though not exclusively by Swift), coupled with a community- 
wide concerted eff'ort to discover and study their afterglows, led 
to a substantial increase in the sample of events over the past 
5 years (e.g.. Figures [2l and O. Studies of this sample have 
led to a re-evaluation of th e host galaxy demographics and the 
redshift distribution (e.g., (JBerger et al.L I2OO7I) ). In particular, 
as of late 2010, the sample of short GRBs with X-ray detections 
(positions of ~ 2 - 5'' radius) numbers about 40. Of these, 
about 20 events have been detected in the optical/UV/near-IR 
and/or radio, leading to positional uncertainties of ~ 0.1 - 0.5''. 
Host galaxies have been identified for nearly all of the bursts 
with sub-arcsecond positions (15/20), and putative hosts have 
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Figure 3: Ground-based images of several short GRB hosts obtained with the 
Magellan 6.5-m telescopes and the Gemini 8-m telescopes. All images are 
20" on a side, with the exception of GRB 060502b which is twice as large. 
The large circles mark the XRT error regions, while smaller circles mark the 
positions of the optical afterglows (when available). Arrows mark the positions 
of the identified hosts. From iBerger et al.i t2007i) . 



also been identified for a substantial fraction of the bursts with 
only X-ray positions (when deep searches have been made). At 
the present, 16 redshifts have been measured between the two 
samples (Figures |4] and O. 

The events with only X-ray positions have two shortcomings. 
First, the probability of chance coincidence for the typical host 
magnitudes within the XRT error circles (~ 21-26 mag) is 
~ 10"^ - 1. Second, in some cases there is disagreement about 
the position and radius of the XRT error circles between various 
groups, leading to systematic uncertainties in host associations. 
Luckily, in the subsequent discussion of detailed host proper- 
ties no substantial diff'erence in the sample with and without 
optical afterglows is found, suggesting that any spurious galaxy 
associations are at most a minor contaminant. 

Using optical follow-up observations of nine short GRBs 
with X-ray and/or optical afterglows (available by the end of 
2006), IBerger et all (I2007h found that eight of the nine bursts 



were likely associated with much fainter galaxies (7? ~ 23 - 26 
mag) than the first few events. By comparison to this early 
sample (with R ~ 17-22 mag and z ^ 0.5), as well as 
the hosts of long GRBs and large field galaxy samples, it was 
demonstrated that these new host galaxies likely reside at z ~ 1 
and beyond. A specific early case for a z ^ 1 origin was 
GRB 060121 based on af t erglow photometric redshi ft estimates 
dde U garte Postigo et al.[ I2OO6I: iLevan et al.1. l2006h . Spectro- 
scopic redshifts for the four brightest galaxies in this expanded 
samp le led to measurements of z ~ 0.4-1.1 (IBerger et al. 
2007k : see Figure IH Subsequent observations have con 




5400 5800 6200 6600 7000 7400 
Wavelength (A) 



[ON] GRB 061210 = 
o 




Call 



GRB 061217 

F([OII])- 1.0x10"^ 
SFR-2.5 



[ON] 



#V%|M^ 



iki- 



5000 5500 6000 6500 7000 
Wavelength (A) 

15 



6000 6200 6400 6600 6800 7000 
Wavelength (A) 



JO 

E 10 

Q) 

00 

7 5 



GRB 060801 

F([OII]) - 1.3x1 0" 
SFR-6.1 



[ON] 



I o^H»%^^ 



m 




6500 7000 7500 8000 

Wavelength (A) 



5500 6000 6500 7000 7500 
Wavelength (A) 



Figure 4: Optical spectra of several short GRB host galaxies. The relevant emis- 
sion lines are marked and lead to redshifts of z ^ 0.4 - 1 . 1 . Also indicated are 
the st ar forma tion rates inferred from the luminosity of the [Oll]A3121 doublet. 
From lBereeJ fmO^ . 
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Figure 5: The redshift distributions of long (gray) and short (red) GRBs as of 
late 2010. The cross-hatched region indicates the redshifts for short GRBs with 
sub-arcsecond positions, while the open histogram includes the redshifts for 
host galaxies identified in some XRT error circle (~ 2 - 5" radius). 
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Figure 6: X-ray luminosity of short GRB afterglows normalized to t = I day (a 
proxy for the afterglow kinetic energy, 6eEK,iso) plotted as a function of £'y4so 
for events with a known redshift (solid black circles), redshift constraints (open 
black circles), and without any redshift information (gray symbols connected 
by dotted lines for z = 0.1 and 1, corresponding roughly to the lowest and high- 
est redshifts securely measured to date). The isotropic-equivalent relativistic 
energies are as least as high as 10^^ erg, and may approach ~ 10^^ erg for some 
short bursts. From iBergen t2007i) . 



redshift distribution (in comparison to that of long GRBs) is 
shown in Figure [3 

One of the crucial ramifications of the measured redshift dis- 
tribution is the energy budget of the y-ray emission and blast- 
wave. In Figure Si show the isotropic-equivalent afterglow X- 
ray luminosities of short GRBs as a function of their isotropic- 
equivalent y-ray energies. Th e forrner is a proxy for the af- 
terglow kinetic ene rgy (Kuman. 120001: iFreedman and WaxmanL 
I2OOII: iBerger et all l2003a). I find that both quantities span 
several orders of magnitude, with £^y,iso ~ 10^^ - 10^^ erg 
(although most short bursts have values of 10^^ - 10^^ erg) . 
This range i s similar to that for long GRB s (IFrail et al .11200 U 
Bergeret aD, [2003a. b ; Bloom et al.l l2003h . and indicates that 



either short GRBs can truly produce a broad range of energies, 
or instead exhibit a wide range of collimation angles. Due to 
the general faintness of short GRB afte rglows, strong evidence 
for collimation exist s in only one case dSoderberg et al.L 2006; 
Burrows et al.Ll2006b . but additional cases are possible (JBergei . 



firmed a broad range of redshifts (e.g ., ( Graham et al.l l2009l : shifts than initially suspected, z ~ 1 (IBerger et al.L l2007b . and 



3. Host Galaxy Luminosities, Metallicities, and Star For- 
mation Rates 

The secure association of at least one short GRB (050724) 
with an elliptical galaxy ("Berger et all 120051: iBarth elmy et al.l 
12005; Gorosabel et al., 2006 ) demonstrated unambiguously that 
some of the progenitors are related to an old stellar popula- 
tion. However, as discussed in the previous section, a substan- 
tial fraction of short GRBs (1/3 - 2 /3) reside at higher red 

kT ■ 



Antonelli et al.l 120091: iLevesque et al.l I2OIOI) ), and the current spectroscopic observations indicate that most of these galaxies 



are undergoing active star formation (Figu re O (IBer ger et al 
20071: ID' Avanzo et all Eo09l: iGraham et all [2009)). Indeed, in 



the sample of short GRBs localized to better than a few arc- 
seconds about 50% of the bursts occur in star forming galaxies 
compared to only ^ 10% in elliptical galaxies; the remaining 
^ 40% are currently unclassified due to their faintness, a lack of 
obvious spectroscopic features, or the absence of deep follow- 
up observations. This result raises the question of whether some 
short GRBs are related to star formation activity rather than to 
an old stellar population, and if so, whether the star formation 
properties are similar to those in long GRB host galaxies. The 
answer will shed light on the diversity of short GRB progeni- 
tors. 
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Figure 7: Specific star formation rates as a function of redshift for the host 
galaxies of short GRBs (black squares), long GRBs (blue circles) and field 
galaxies from the GOODS-N survey (red crosses; (Kobulnicky and Kewley, 
[2004)). Upper limits for the elliptical hosts of GRBs 050509B and 050724 are 
also shown. The cross-hatched region marks the median and standard deviation 
for the long GRB host sample. The inset shows the cumulative distributions for 
the three samples. The K-S probability that the short and long GRB hosts are 
drawn from the same distribution is only 0.3%, while th e strong overl ap with 
the field sample leads to a K-S probability of 60%. From lBergeJ ( l2009[) . 



In the following discussion, I compare several aspect of short 
and long GRB host galaxies: Luminosities, metallicities, and 
star formation rates. A comparison of the masses and stel- 
lar population ages is carried out in the subsequent section. 
For the current sample of short GRB hosts, the distribution 
of absolute rest-frame 5-band magnitudes {Mb) rang e s from 
about 0.1 to a few L, (Berger, 2009; Prochaska et al.1 l2006h . 
The star formation rates (mostly infer red from the [011 1/13727 
line using the standard conversion dKennicuttl [l998j; Fig- 
ure 131) ran ge froni about 0.1 to 10 Mq yr~^ for the star form- 
ing hosts (IBergen. l2009l: Id' Avanzo et al.L l2009b . In the case 
of the elliptical hosts the uppe r limits are < 0.1 Mq yr~ ' 



iBerger et al.1 l2005l:lBloom et al.1 [2006: Prochaska et al.ll2QQ6l : 



iBergeii l2009b . Combined with the absolute magnitudes, the 
specific star formation rates (SSFR) are SFR/Lg ^1-10 M© 
yr"^ L~^ for the star forming hosts, and < 0.03 M© yr"^ L~^ for 
the elliptical hosts. The SSFR values as a function of redshift 
are shown in Figure [TJ 



For fivqj host galaxies in the current sample, there is also 
sufficient spectral information to measure the metallicity 
Berge^ (I2Q09I) : IProchaska elaP (2006); D' Avanzo et al. 
(■2009b . I use the standard metallicity diagnostics . 



/^23 = [(-^rOIIU3727 + -^rQI I IU^495 9.5007)/-^H/3l dPagcl Ct al 



1979|; [Kobulnicky and Kewleyl 120041) and i^[Nii]/i6584/^Ha. The 
value of /?23 depends on both the metallicity and ionization 
state of the gas, which is determined using the ratio of oxygen 
lines, O32 = /^[Oiii]i^4959,5007/^[oii]i3727. I notc that the R23 
diagnostic is d ouble- valued with l ow an d high metallicity 
branches (e.g., (IKewley and DopitaL l2002b ). This degeneracy 
can be broken using the ratio [NIIJ/Hof when these lines are 
accessible. To facilitate a subsequent comparison with field 
ga laxy samples I use the 7? 2 3, 0?>2 , and [NII]/Hck calibrations 
of iKobulnicky and Kewleyl (|2004|) . The typical uncertainty 
inherent in the calibrations is about 0.1 dex. 



Adopting the solar metallicity fr om lAsplund et al.l (l2005h 



12 -h log(0/H) = 8.66'Berger| (l2009h find 12 -h log(0/H) ^ 8.6 
for the upper R22, branch and ^ 8.0 - 8.5 for the lower branch 
for the host of GRB 061006. For the host of GRB 070724 they 
find 12 -h log(0/H) ^ 8.9 for the upper branch, and ^ 7.6 - 8.1 
for the lower branch. A similar range of values is found for 
the host of GRB 061210, but the ratio F[Nii]//^Ha ~ 0.2, indi- 
cates 12 -h log(0/H) > 8.6, thereby breaking the degeneracy 
and leading to the upper branc h solution, 12 -h log(Q/H ) ^8.9. 
For the host of GRB 05 1221 A dSoderberg et all (20061) . similar 
values to those for the host of GRB 070724 are inferred. Fi- 
nally, for the host galaxy of GRB 050709, the [Nil] /Ha ratio 
indicates 12 -h log(0/H) ^ 8.5. The dominant source of un- 
certainty in this measurement is the unknown value of O32, but 
using a spread of a full order of magnitude results in a metal- 
licity uncertainty of 0.2 dex. For the hosts with double- valued 
metallicities (GRBs 05 1221 A, 061006, and 070724) I follow 
the conclusion for field galaxies of similar luminosities and red- 
shifts t hat the appropriate values are t hose for the R22, upper 
branch ("Kobulnickv and Kewlev', 2004 ). Thi s conclusions was 
advocated by Kobulnicky and Kewleyl (l2004l) based on galaxies 
in their sample with measurements of both 7?23 and [Nil] /Ha. 
It is similarly supported by our inference for the host galaxy of 
GRB 061210. Future near-IR spectroscopy covering the [Nil] 
and Ha lines will test this hypothesis. The metallicities as a 
function of host luminosity are shown in Figure (8] 

To place the host galaxies of short GRBs in a broader con- 
text I compare their properties with those of long GRB hosts 
and field star form ing galaxies from the GOODS-N survey 
dKobulnicky and K ewley. 2004). In terms of absolute magni- 
tudes, the long GRB hosts range from Mb ~ -16 to -22 mag, 
with a media n value of {Mb) ~ -19.2 mag {{Lb) ~ 0.2 L*; 
JBerger et al.[ i2007) ). Thus, the long GRB hosts extend to 
lower luminosities than the short GRB hosts, with a median 
value that is about 1.1 mag fainter. A K-S test indicates that the 
probability that the short and long GRB hosts are drawn from 



^A sixth host, GRB 071227, has an inferred metallicity of about 0.4 - 1 
Z0, but this was inferred in the absence of detected hydrogen B aimer lines 
JD' Avanzo et all [2009) . and the values are therefore prone to large systematic 
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Figure 8: Metallicity as a function of 5-band absolute magnitude for the host 
galaxies of short GRBs (black squares) and long GRBs (blue circles). The gray 
bars mark the 14 - 8 6 percentile ran ge for galaxies at z ~ 0.1 from the Sloan 
Digital Sky Survey jTremonti et al.L |2004), while red crosses designat e the 
same field galaxies at z ~ 0.3 - 1 shown in Figure 0(Kobulnickv and KewlevI, 
[2004). Both field samples exhibit a clear luminosity-metallicity relation. The 
long GRB hosts tend to exhibit lower than expected metallicities ( Stanek et al., 
|2006), while the hosts of short GRBs have higher metallicities by about 0.6 
dex, are moreover in excellent agreement with the luminosity-metallicity rela- 
tion. From Bergej i2009l) . 



the same underlying distribution is 0.1. On the other hand, a 
comparison to the GOODS-N sample reveals a similar distri- 
bution, and the K-S probability th at the short G RB hosts are 
drawn from the field sample is 0.6 (lBergeill2009b . 

A similar conclusion is rea ched based on a comparison of 
specific star formation rates (IBergen. l2009h . For long GRB 
hosts the inferred star formation rates range from about 0.2 to 
50 Mo yr"\ and their specific star formation rates are about 
3 - 40 Mq yr~^ L~^, with a m edian value of about 10 M© yr"^ 
L~^ (IChristensen et al.L l2004l) . As shown in Figure [71 the spe- 



cific star formation rates of short GRB hosts are systematically 
below those of long GRB hosts, with a median value that is 
nearly an order of magnitude lower. Indeed, the K-S probabil- 
ity that the short and long GRB hosts a re drawn frorn the same 
underlying distribution is only 0.003 (JBergeii l2009h . This is 
clearly seen from the cumulative distributions of specific star 
formation rates for each sample (inset of Figure [7]). On the 
other hand, a comparison to the specific star formation rates of 
the GOODS-N field galaxies reveals excellent agreement (Fig- 
ure [7]). The K-S probability that the short GRB hosts are drawn 
from the field galaxy distribution is 0.6. Thus, short GRB hosts 
appear to be drawn from the normal population of star forming 
galaxies at z ^ 1, in contrast to long GRB hosts, which have 
elevated specific star formation rates, li kely as a result of pref- 
erentially young starj burst populations JChristensen et al.ll2004l : 



Savaglioetalll2008h . 



Finally, the metallicities measured for short GRB hosts are in 
excellent agreement with the luminosi ty-metallicity relation for 
field galaxies at z ~ 0.1 - 1 (Figure[8l (JKobulnicky and KewleyL 



2004 iTremonti et al.l. 120041) ). The two hosts with Mg ^ -18 
mag have 12 + log(0/H) ^ 8.6, while those with Mb ~ -20 to 
-21 mag have 12 + log(0/H) ^ 8.8 - 8.9, following the general 
trend. On the other hand, the short GRB host metallicities are 
systematically higher than those of long GRB host s, which have 



been argued to have lower than expected values (Stane k et al 



20061). The median metallicity of short GRB hosts is about 0.6 
dex higher than for long GRB hosts, and t here is essenti ally no 
overlap between the two host populations (lBergeiil2009h . 

To conclude, the short GRB host sample is dominated by 
star forming galaxies, but these galaxies have higher luminosi- 
ties, lower star formation rates and specific star formation rates, 
and higher metallicities than the star forming host galaxies of 
long GRBs. Instead, the short GRB host sample appears to be 
drawn from the field galaxy population. These results suggest 
that while short GRB hosts are mainly star forming galaxies, the 
progenitors most likely trace stellar mass rather than the modest 
on-going star formation activity. 

4. Host Galaxy Stellar Masses and Ages 

To more comprehensively address whether short GRBs trace 
stellar mass alone (as would be expected for an old progeni- 
tor population), it is essential to determine the stellar masses 
and population ages of short GRB host galaxies, primarily in 
comparison to the general galaxy st ellar mass function This 
analysis was recently carried out by Leibler an d Berger (2010) 
using multi-band optical and near-IR data for 19 short GRB 
hos ts. The resulting spectral energy distributions were fit with 
the iMarastonI (|2005) stellar population models to extract two 
crucial parameters: stellar mass (M^) and population age (r). 
The range of possible masses was assessed using three ap- 
proaches. First, using single stellar population (SSP) fits, which 
provide an adequate representation for the early-type hosts, but 
tend to under-estimate the total mass and population age of star- 
forming hosts. At the other extreme, the near-IR data alone 
were modeled with a stellar population matched to the age of 
the universe at each host redshift. This approach uses the maxi- 
mum possible mass-to-light ratio to extract a maximal mass for 
each host galaxy. Finally, as a more realistic estimate for the 
star-forming hosts, hybrid young+old stellar populations were 
used. Examples of all three approaches are shown for the host 
of GRB 050709 in Figure[9l 

The resulting mass distributions are shown in Figure [TOl For 
comparison I also present the mass distributions for long GRB 
hosts, which were analyzed with the same models for the pur- 
pose of a uniform comparison. The SSP masses span three or- 
ders of magnitude, Mssp ~ 6 x 10^ - 4 x 10^^ M©, with a 
median value of (Mssp) ~ 1.3 x 10^^ M©. Dividing the sam- 
ple into early- and late-type host galaxies, the former span the 
range Mssp ~ (2 - 40) x 10^^ M©, while the latter have much 
lower masses of Mssp ~ (0.06 - 2) x 10^^ M©. The clear dis- 
tinction between the two samples partially reflects the bias of 
single age SSP models, which for the late-type hosts are domi- 
nated by the young stellar population and hence under-estimate 
the contribution of any older stellar populations. For the max- 
imal masses the range is Mjviax ~6xlO^-8xlO^^ M©. The 
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Figure 9: Optical and near-IR spectral energy distribution of GRB 050709 with the three models used to extract the stellar mass and population age. Data are shown 
as black circles with error bars, and synthesized model fluxes are shown as red squares. Left: Single age SSP model. Center: Young+Old SSP model (magenta=old; 
cyan=young) with the old population age fixed at the age of the universe at the redshift of the burst (z = 0.161 in this case). Right: Maximal mass model with the 
population age fixed at the age of the universe and using only the /T-band photometry. The Young+Old model leads to total masses intermediate between the single 
age SSP and the maximal models, and has younger ages for the young population than the single age SSP model. From Leibler and Bergeil ( 120 IQl) . 



median mass is (Mmex) ~ 1x10^^ M©, about an order of magni- 
tude larger than for the single age SSP masses, and only slightly 
larger than the stellar mass of the Milky Way. As expected, the 
ratios of maximal to SSP masses for the early-type hosts are 
modest, MMax/^ssp ~ 2 - 8, since these hosts are already dom- 
inated by old stellar populations. However, for the late-type 
hosts the corrections are significant, Mmex/^ssp ~ 5 - 60, with 
a median ratio of about an order of magnitude. Finally, using 
the young-hold models for the late-type hosts, and the single 
age SSP values for the early-type hosts, the resulting masses 
are M ^ 2x10^-4x10^^ M©, with amedian of <M) ^ 5x10^^ 
Mo. 

Stellar population ages can only be inferred for the single 
age SSP models since for the maximal and young-hold models 
the inherent assumption is a population with the age of the uni- 
verse at each host redshift. The distribution of ages is shown 
in Figure [TT] with the values ranging from about 30 Myr to 4.4 
Gyr. The median age is (tssp) ~ 0.3 Gyr for the full sample, 
with (tssp) ~ 0.25 Gyr for the subset of late-type hosts and 
('Tssp) ~ 3 Gyr for the subset of early-type hosts. 

The long GRB hosts, on the other hand, have lower masses 
and younger population ages. For the SSP model the mass 
range is Mssp ~6xlO^-2xlO^^Mo, with a median value 
of (Mssp) ~ 1.3 x 10^ M© (Figure [T0|. The maximal masses 
span MMax ~ 9 X 10^ - 9 X 10^^ M©, with a median value of 
(MMax) ~ 4.0 X 10^ Mo (Figure [T0|. The SSP stellar popu- 
lation ages span about 10 to 570 Myr, with a median value of 
(tssp) ~ 65 Myr. 

4.1. Host Demographics 

In the redshift range relevant for short GRBs (z ~ 0.2 - 1) 
roughly an equal fraction of the cosmic ste llar mass budget r e- 



sides in early- and late-type galaxies (e.g., (Illbert et al.Ll2010l) ) 



Therefore, if short GRBs track stellar mass alone we expect 
a roughly one-to-one ratio of galaxy types. This does not ap- 
pear to be the case. For example, within the sample of short 
GRBs with optical afterglows (20 events), only 2 are unam- 
biguousl y hosted by early-t ype galaxies (GRBs 050724 and 
100117; (JBerger et al.l l2005h . Fong et al. in prep.), while 8 are 
unambiguously hosted by late-type galaxies; the probability of 



obtaining this ratio from an intrinsic one-to-one distribution is 
only 0.04. The identity of the remaining 9 hosts is unclear at the 
present due to their faintness or the lack of underlying galaxies 
at the burst positions. Still, unless nearly all of these bursts 
were hosted by early-type galaxies, the resulting ratio appears 
to be skewed in favor of late-type host galaxies with on-going 
star formation activity. The same result holds true even when 
considering the bursts with only X-ray afterglow positions and 
identified hosts. 

Thus, the host galaxy demographics suggest that short GRBs 
do not track stellar mass alone, or phrased alternatively, they 
do not have a delay time distribution that is skewed to old ages 
of ~ few Gyr. It is possible, however, that this result is in- 
fluenced by secondary factors such as the typical circumburst 
density or intrinsic difl'erences in the energy scale and after- 
glow brightness as a function of galaxy type (possibly remi- 
niscent of the difl'erences in peak lumino sity for Type la su- 
pernovae in early- and late-type galaxies; ('Hamuy et al.l l2000c 
■Mannucci et al.ll2006b ). If such difl'erences lead to fainter after- 
glows (or prompt emission) for short GRBs in early-type galax- 
ies, this would suppress the early-type fraction. Although the 
modest size of the host sample, and the substantial fraction of 
short GRBs with only y-ray positions (~ 1/3 of all events), pre- 
vent definitive conclusions, it does not appear that the optical 
afterglo ws of short G RBs in early- and late-type galaxies are 
distinct dBergeilEoTQl) : see ^ 



4.2. Comparison to the Galaxy Mass Function 

I next turn to a compariso n of the inferred stellar masses 
with the galaxy mass function (ILeible r and Berger, 2010). The 



cumulative distribution of stellar masses for the short GRB 
hosts is shown in Figure [12] with the range of possible masses 
bounded by the single age SSP and maximal values. I also 
present a breakdown of the sample into early- and late-type 
galaxies, each spanning the same range. For the late-type hosts, 
the intermediate young-hold values are also shown. To com- 
pare these distributions to the distribution of galaxy masses I 
also plot the cumulative stellar mass function weighted by mass, 
i.e., the fraction of stellar mass in galaxies above some mass. 
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Figure 10: Top: Histograms of inferred stellar masses from single stellar popu- 
lation fits for the hosts of short (black) and long (gray) GRBs. The inset shows 
the cumulative distributions, including for the subset of late-type short GRB 
hosts (blue). The median values for the three samples are given in parentheses, 
and the Kolmogorov-Smirnov probabilities that the distributions of short and 
long GRB hosts, as well as star forming short GRB and long GRB hosts are 
drawn from the same distribut ion are provided in t he ins et. Bottom: Same but 
for the maximal masses. From lLeibler and Bergej ( |201Q|) . 



f(> M), given by the equation: 

f(>M)=^ 

where 0(M) is the Schechter mass function: 



M 



0(M) = 0*- exp--). 
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Several determin ations of 0(M) a re used for this compar- 
ison, including the ICole et alJ (1200 I j mass function from the 
2MASS/2dF catalogs for all galaxy ty pes at z - (M* = 10^^-^^ 

m 



Short GRBs: SSP (-0.5) 
Short GRBs: late (-0.6) 
Long GRBs: SSP (-1.2) 




9 
8 

7 

6 

^_ 

CD 

-Q 5 
Z 4 



-1 -0.5 
log [age / Gyr] 



Figure 1 1 : Histograms of inferred stellar population ages from single stellar 
population fits for the hosts of short (black) and long (gray) GRBs. The in- 
set shows the cumulative distributions, including for the subset of late-type 
short GRB hosts (blue). The median values for the three samples are given in 
parentheses, and the Kolmogorov-Smirnov probabilities that the distributions 
of short and long GRB hosts, as well as star forming short GRB and long GRB 
hosts are drawn from the same distribution are provided in the inset. From 
iLeibler and Bereej ( l20icil) . 



Mo, a = -1.16); the 'Bell et all (l20Q3l ) mass function for 
late-type galaxies from 2MASS/SDSS converted to a Salpeter 
IMF for comparison with our inferred value s (M* = 10^^-^^ 
Mq, a = -1.27); and the Ilbert et al.l (120101) mass functions 
from the COSMOS survey for quiescent galaxies at z ~ 0.3 
(M* = 10^^-^^ Mo, a = -0.91) and intermediate- activity galax- 
ies at z ~ 0.5 (M* = lO^O-^^ Mo, a = -1.02), matched to the 
redshifts of the early- and late-type short GRB host galaxies in 
our sample. The resulting distributions of /(> M) for the vari- 
ous mass functions are shown in Figure [T2l 

The agreement (or lack thereof) between the short GRB 
hosts and the galaxy mass functions is assessed using the 
Kolmogorov-Smirnov (K-S) test. For the full sample there is 
negligible probability that the distribution of single age SSP 
masses is drawn from the galaxy mass function, with P ^ 
8 X 10"^. On the other hand, for the maximal mass distribution 
the probability is P ^ 0.6 indicating that for these masses the 
short GRB sample is fully consistent with the galaxy mass func- 
tion. Using the intermediate case of SSP masses for the early- 
type hosts and the young -hold masses for the late-type hosts, the 
probability is P ^ 0.3, indicating that this combination is also 
consistent with the galaxy mass function. 

Separating th e early-type hosts , their SSP masses are con- 
sistent with the lllbert et al.l (120101) mass function of quiescent 
galaxies (P ^ 0.8); their large maximal masses, on the other 
hand, are inconsistent with the mass function, with P ^ 0.007. 
Finally, for the late-type hosts there i s a clear inconsiste ncy of 
the single age SSP masses with the lllbert et al.l (12010 1) mass 



Mo, a = -1.18); the nearly identical iPanter et al.l (120041) mass 
function from SDSS for all galaxy types at z ~ (M* = 10^^-^^ 



function of intermediate- activity galaxies, with P ^ 4x 10"^. 
However, the maximal mass distribution is consistent with the 
mass function (P ^ 0.3), while the young -hold mass distribution 
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Figure 12: Cumulative distributions for the full sample of single age simple 
stellar population (SSP) masses, maximal masses, and combined young+old 
and SSP masses for the late- and early-type hosts, respectively (black; bottom 
panel). The upper panel shows a breakdown by galaxy type (late-type: blue; 
early-type: red). The shaded regions represent the range of possible stellar 
masses since the SSP masses, which are effectively light-weighted values, are 
most likely an under-estimate, while the maximal masses make the extreme 
assumption that all hosts are dominated by populations with the age of the uni- 
verse. For the late-type hosts the total masses from a young-i-old SSP fit are 
also shown; these are more closely representative of the total mass. Also shown 
are the fractions of total stellar mass in galaxies with mass, > M, calculated 
from several published galaxy stellar mass functi ons at z ~ - 2 (cyan and ma- 
genta lines; (Cole et al., 2001; Bell et al. j2003l : lPanter et al.L[20Q4l : lllbert et SI. 
I2OIO)); for the Ilbert et al. (2010) mass function the z ~ 0.5 bin is used, ap- 
propriate for the short GRB sample, and separately plot the mass function for 
quiescent galaxies and for intermediate-activity galaxies, which resemb le the 
intermediate star formation activity in short GRB hosts dBergej . 12009). The 
comparison indicates that short GRBs trace galaxy mass only if the bulk of the 
late-type hosts have close to maximal masses. The subset of early-type hosts 
appears to fai thfully trace the mass fun ction of galaxies for the SSP-derived 
masses. From lLeibler and Bergej ( l201Qh . 



is marginally consistent {P ^ 0.1). 

To summarize, the distribution of short GRB host masses is 
compatible with the overall mass distribution of galaxies only 
if their stellar masses are given by the SSP masses for the early- 
type hosts and the maximal or young+old masses for the late- 
type hosts. Since the opposite scenario (maximal masses for 
the early-type hosts and SSP masses for the late-type hosts) 
is unlikely, the existing sample of short GRB hosts is consis- 
tent with the galaxy mass function (ILeibler and BergenbOlQI) . 
Equivalently, this means that short GRBs may indeed track stel- 
lar mass alone. However, I caution that the host demographics 
seem to be at odds with the expected equal fractions of total 
stellar mass in early- and late-type galaxies, unless nearly all 



of the unidentified hosts are early-type galaxies. This, along 
with the somewhat lower than expected masses of the late-type 
hosts, leaves open the possibility that at least a subset of short 
GRB progenitors track star formation activity rather than stellar 
mass. 

4.3. Comparison to Long GRB Hosts 

Despite the possibility that some short GRB progenitors may 
track star formation activity, the inferred stellar masses and 
population ages of short GRB hosts are generally distinct from 
those of lon g GRB hosts in both the si ngle age SSP and maxi- 
mal models (ILeibler and Bergen bOld) . Most importantly, this 
is true for the subset of late-type hosts. In the framework of the 
single age SSP model the K-S probability is only 0.006 that the 
long and short GRB hosts are drawn from the same mass dis- 
tribution. The probability is higher for the subset of late-type 
short GRB hosts, P ^ 0.1 (Figure[TO|. However, since the SSP 
values represent the mass of only the young stellar populations, 
they are mostly reflective of the star formation activity rather 
than the total stellar mass. Using instead the maximal masses, 
the K-S probability that the long GRB hosts and late-type hosts 
of short GRBs are drawn from the same sample is negligible, 
P ^ 4x 10"^ (Figure [T0|, demonstrating that they are distinct 
galaxy populations. A similar conclusion is apparent from a 
comparison of the single age SSP population ages. The K-S 
probability that the long GRB hosts and late-type hosts of short 
GRBs are drawn from the same distribution is only P ^ 0.006. 

Thus, the long GRB hosts have significantly lower stellar 
masses than the subset of late-type short GRB hosts, and their 
young stellar population are significantly younger. Indeed, 
a comparison of the long GRB host maximal masses to the 
Ilbert et al.l (120101) mass function of high-activity galaxies at 
z ~ 0.7 (appropriate for the long GRB sample considered here) 
indicates a K-S probability of only 0.002 that the long GRB 
hosts are drawn from the galaxy mass function. This is consis- 
tent with our understanding that their massive star progenitors 
select galaxies by star formation (and perhaps additional factors 
such as metallicity), but not by stellar mass. 

The apparent distinction between long GRB hosts and the 
late-type hosts of short GRBs in terms of their stellar masses 
and young population ages strengthens the conclusion in O 
based on the star formation ra tes, specific st ar formation rates, 
luminosities, and metallicities (Berger, 2009; . In essentially ev- 
ery property the late-type short GRB hosts point to a population 
of more quiescent, massive, and evolved galaxies than the hosts 
of long GRBs. I therefore conclude that this rules out the idea 
that short GRB progenitors in late- type hosts are mass ive stars 
identical to long GRB progenitors (Ivirgili et al.L l2009k . even if 
the short GRBs in late-type galaxies indeed track star formation 
rather than stellar mass. 

4.4. The Delay Time Distribution 

A determination of the delay function from the derived stel- 
lar population ages is complicated by two primary factors. First, 
they rely on the assumption that the short GRB progenitors in 
each host were formed within the inferred stellar population. 



This assumption is justified statistically both for an association 
of the progenitors with stellar mass and with star formation ac- 
tivity, as long as we can appropriately normalize the rates of 
short GRBs. Second, while it is possible to determine single age 
SSP ages from the broad-band photometry, these data are not 
sufliicient to provide an age breakdown (by mass) for multiple 
stellar components. Indeed, for the hybrid young-Fold model, 
the age of the old populati on has to be fixed (to the age of the 
universe in the analysis of Leib ler and Bergej (l201Qh ). Still, in 
the young-Fold model, the bulk of the mass (^ 55 - 99%) is 
contained in the old stellar population, and so the progenitors 
would have "old" ages (r > tssp) if they tracked stellar mass. 

As a result of these limitations it is only possible to explore 
the implications of two main scenarios, namely that short GRBs 
track mass and/or star formation activity. In the context of the 
former scenario I have shown in ^4.21 that the short GRBs in 
early-type hosts trace stellar mass. Therefore, their SSP ages 
can be used to provide a rough estimate of the progenitor ages, 
T ^ 0.8 - 4.4 Gyr, with a median of about 3 Gyr. On the other 
hand, for the late-type hosts (for which no credible information 
on the mass- weighted stellar population age can be extracted), it 
is possible to infer a typical delay relative to the most recent star 
formation episode under the assumption that these progenitors 
track star formation activity. I find SSP ages of r ^ 0.03 - 0.5 
Gyr, with a median of about 0.25 Gyr, or young-hold ages of 
about 0.01 - 0.4 Gyr with a median of about 0.1 Gyr. 

Thus, if short GRBs follow both stellar mass (in early-type 
galaxies) and star formation activity (in late-type galaxies), the 
typical delay times are about 3 and 0.2 Gyr, respectively. 

5. Offsets and the Sub-Galactic Environments 

I next turn from a galactic- scale investigation of short GRB 
host environments to the sub-galactic scale. In general, the 
sub-galactic environments of cosmic explosions provide pow- 
erful insight into the nature of their progenitors. For exam- 
ple, in the case of long GRBs, the distribution of projected 
physical and host-normalized offsets relative to the host centers 
matched the expected distribution for massive stars in an ex- 



ponential disk (Bloom et al.U2002b . Moreover, long GRBs are 



spatially correl a ted wi th the brightest UV regions of their hosts 
(JFruchter et al.', zOOo). Both of these studies have relied on 
high angular resolution Hubble Space Telescope (HST) obser- 
vations to spatially resolve the hosts and astrometrically locate 
the GRB positions to pixel-scale accuracy. As a by-product, 
these observations also provided detailed morphological infor- 
mation for the hosts (e.g ., Sersic index, eff'ective radius; e.g. 
dWainwright et al.','2007^). 

Individual off^sets have been measured for several shor t 
GRBs (e.g., dBerger et al.l. l2005l: IFox etaP. l2005l: iBloom et al 



20061) ), and a study relying on ground based dat a without a 
compl ete astrometric treatment was published by iTroja et aL 



( 2008k . However, the first systematic study of short GRB off- 
sets, as well as their galactic environments and host morpholo- 
gies was recently published by iFong et al.l ( [2010 ). The sam- 
ple includes ten short GRBs (spanning May 2005 to December 
2006), of which seven have been localized to sub-arcsecond 



precision, and of those, s ix are robustly as sociated with host 
galaxies (for details see (JFong et al.l l2010h ). Illustrative ex- 
amples of HST host images and morphological model fits are 
shown in Figure [131 
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Figure 13: Top-left: //577WFPC2/F814W image of the host galaxy of 
GRB 05 1221 A with a 5cr error circle representing the afterglow position. Top- 
center: Sersic model fit from galf it. Top-right: Residual image. Bottom: 
Same, but for the host galaxy of GRB 061006. From iFong et al.i (20101) . 



5.7. Morphological Analysis 

The two-dimensional surface brightness profiles were used 
to determine the hosts' eff'ective radii and morphological prop- 
erties such as the Sersic n index. These quantities are crucial 
for a comparison of the morphologies and sizes to those of long 
GRBs, as well as for normalization of the projected off'sets rel- 
ative to the galaxy size. Three hosts (GRBs 050709, 05 1221 A, 
and 061006) are best modeled with n ^ I, corresponding to an 
exponential disk profile, while two hosts (GRBs 050509B and 
050724) are best modeled with n ^ 3 and ^ 5.6, respectively, 
typical of elliptical galaxies. These classifications are in perfect 
agreement with their spectroscopic properties. The final three 
hosts in the sample (GRBs 051210, 060121, and 060313) are 
faint, and as a results can be modeled with a wide range of n 
values, although ^ ~ 1 is preferred in all three cases. Therefore, 
of the eight short GRB hosts only two can be robustly classi- 
fied as elliptical galaxies based on their morphology. A similar 
fraction was de termined inde pendently from spectroscopic ob- 
servations (gl (!Berg ei^j2009)). 

The morphological analysis also yields values of the galaxy 
eff'ective radii, r^. A range of ^ 0.2 - 5.8'' is found, correspond- 
ing to physical scaleq^ of about 1.4 - 21 kpc. The smallest 
eff'ective radius is measured for the host of GRB 060313, while 
the elliptical host of GRB 050509B has the largest eff'ective ra- 
dius. The median value for the sample is r^ ^ 3.5 kpc. The 
eff'ective radii as a function of n are shown in Figure \T5\ Also 



^For the faint hosts without a known redshift (GR Bs 051210, 060121 , 
060313, and possibly 061201) it is assumed that z = 1 ^Berger et all l2007h . 
and take advantage of the relative flatness of the angular diameter distance as a 
function of redshift beyond z ~ 0.5. 
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Figure 14: One-dimensional radial surface brightness profiles for short GRB 
host galaxies derived using the IRAF task ellip se. The gray lines are Sersic 
model fits to the surface brightness profiles. From lFong et al.l J2010l) . 
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a similar analysis carried out by lWainwright et alJ (12007^ . 

Two clear trends emerge from the morphological comparison 
of short and long GRB hosts. First, all long GRB hosts have 
n < 2.5, and the median value for the population is (^) ^ 1.1 
( Wainwright et al.L l2QQ7h . Thus, they are all morphologically 
classified as exponential disk galaxies, while 2 of the 8 short 
GRB hosts exhibit de Vaucouleurs elliptical galaxy profiles. 
However, for the short GRB hosts with ^ < 2, the distribu- 
tion of n values ap pears to be similar to that of long GRB hosts 



( Fongetal.l.l2QlQh 



Second, the short GRB host galaxies have larger eff'ective 
radii, with ( re) ^ 3.5 kpc, com pared to (r^) ^1.7 kpc for long 
GRB hosts dWainwright et al.ll2()()7l) . A Kolmogorov-Smirnov 
(K-S) test indicates that the probability that the short and long 
GRB hosts are drawn from the same underlying distribution 
of host galaxy effective radii is only 0.04. Thus, short GRB 
host galaxies are systematically larger than long GRB hosts. 
The larger sizes are expected in the c ontext of the we ll-known 
galaxy size-luminosity relation (e.g., (lFreemanLll9 70)) and the 
higher luminosity of short GRB hosts (gl dBergeri. 120091) ). 

An additional striking difference between the hosts of long 
and short GRBs is the apparent dearth of interacting or irreg- 
ular galaxies in the short GRB sample. Of the eight short 
GRB host galaxies with HST observations only one exhibits 
an irregular morphology (GRB 050709) and none appear to be 
undergoing mergers. In contrast, the fraction of long GRB 
hosts with an irregular or merger/ interaction morphology is 
~ 30-60% dWainwright et"alll2007h . The interpretation for the 
high merger/interaction fraction in the long GRB sample is that 
such galaxies represent sites of intense star formation activity 
triggered by the merger/interaction process, and are therefore 
suitab le for the production of massive stars (Wainwrig ht et al. . 
20071). The lack of morphological merger signatures in the 



Figure 15: Eff'ective radii for the short GRB hosts observed with HST as a 
function of their Sersic n values (Figure QJJ. Also shown are t he data for long 
GRB hosts based on HST observations from the sample of Wainwright et al.l 
§007). The hosts of short GRBs 050509B and 050724 have n values typical of 
elliptical galaxies, but the remaining hosts have a similar distribution to that of 
long GRBs (i.e., a median ofn ~ 1, or an exponential disk profile). On the other 
hand, the hosts of short GRBs are larger by about a factor of 2 than the hosts 
of long GRBs, in agreement with their higher luminosities. From Fong et^l 
(120100 . 



the merger and the occurrenc e of a short GRB is > 1 Gyr (e.g., 
( Barnes and HernquistLll992b ). 



5.2. The Offset Distribution 

The location of each short GRB relative to its host galaxy 
center and its overall light distribution was determined through 
differential astrom etry using optical and near-IR images of 
the afterglowqj (F ong et al.L 120 lOi) . With the exception of 



GRB 050709, whose afterglow is directly detected in HST 
observations, ground-based afterglow images from Magellan, 
Gemini, and the VLT were used. The resulting positional un- 
certainties include contributions from the ground-based to HST 
astrometric tie ((Tgb^hst ~ 10-30 mas), the positional un- 
certainty of the afterglow ((T^,grb ~ 1 - 40 mas for optical af- 
terglows and ^ 1.7-5.8 arcsec for X-ray afterglows), and the 
uncertainty in the centroid of the host galaxy (cre,gai ~ 1 - 20 
mas). The resulting combined offset uncertainties for the short 
GRBs with optical afterglows are < 60 mas, corresponding to 
physical offset uncertainties of < 0.5 kpc; the best-measured 
offsets have uncertainties at the level of tens of pc. These off- 
sets also correspond to about 1 HST pixeL 

Based on the resulting astrometric ties Fong et al. (2010) find 
that the projected offsets of short GRBs relative to their host 
centers range from about 0.12 to 17.7''. The corresponding 
projected physical offsets are about 1-64 kpc, with a median 



short GRB sample indicates that if any of the hosts have un- 
dergone significant mergers in the past, the delay time between 



^Optical afterglows have not been detected for GRBs 050509B and 051210. 
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value of about 5 kpc. The largest offsets are measured for GRBs 
050509B and 051210, but those are based on Swift/XRT posi- 
tions only, with statistical uncertainties of 12 and 18 kpc, re- 
spectively, considering only the bursts with sub-arcsecond af- 
terglow positions, the largest offset is 3.7 kpc (GRB 050709), 
and the median offset for the 6 bursts is 2.2 kpc. In the case of 
GRB 061201 the host association remains ambiguous, but even 
for the nearest detected galaxy the offset is about 14.2 kpc. 




1 10 

Offset (kpc) 

Figure 16: Projected physical offsets for short GRBs (black) and long GRBs 
(gray; (Bloom etal., 2002)). The top panel shows a cumulative distribution, 
while the bottom panel shows the differential distribution taking into account 
the non-Gaussian errors on the offsets. The arrows in the bottom panel mark 
the median value for each distribution. The median value for short GRBs, ^ 
5 kpc, is about a factor of 5 times larger than for long GRBs. The arrows 
in the top panel mark the most robust constraints on the offset distribution, 
taking into account the fraction of short GRBs with only y-ray positions, as well 
as short GRBs for which hosts have been identified within XRT error circles 
(thereby providing a typical range of ~ - 30 kpc). Also shown in the top panel 
are predicted offset distributions for NS-NS binary mergers in Milky Way type 
galaxies based on population synthesis models. Good agreement between the 
observed distribution and models, as well as between the robust constraints and 
models is found. From Fong et al. (2010). 

To investigate the offset distribution in greater detail the 
HST sample was supplemented with offsets for GRBs 070724, 
07122 7, 080905 A and 9 0510 from ground-ba sed observa- 
tions (iBerger et al.l 120091: iRau et al.l 12009; Rowlinsonet al. . 
20101) . These bursts have accurate positions from optical af- 



terglow detections. In the case of GRBs 070724 and 071227 
the optical afterglows co incide with the disks of edge-on spi- 
ral g alaxies (Figure [2 (IBerger et al.[ 120091: IP'Avanzoet al. . 
[2009) ). The offsets of the three bursts are 4.8, 14.8, 18.5, and 
5.5 kpc, respectively. 

There are 7 additional events with optical afterglow identifi- 
cations. Of these, two bursts (070707 and 0707 14B) coincide 



with galaxies dPiranomonte et al.l 120081: iGraham et aP. l2009h . 
but their offsets have not been reported by the authors. Based 
on the claimed coincidence a conservative estimate is < 0.5'', 
corresponding^ to < 4 kpc. For GRB 090426 an offset relative 
to one of the knots in the host galaxy complex was reported 
(ILevesque et al.L |2010), but not relative to the host center. Fi- 
nally, four bursts (070809, 080503, 090305, 090515) do not 
have coincident host galaxies to deep limits; these bursts are 
discussed in detail in ^ 

In addition to the bursts with sub-arcsecond positions, several 
hosts have been identified within XRT error circles in follow- 
up observations (GRBs 060801, 061210, 6 1217, 070429B , 



070729, and 080123; dBerger et al.l 120071: iBergeJ. l2009k ). 



Since the putative hosts are located within the error circles, the 
offset s are consistent with zero or may be as large as ~ 30 kpc 
(e.g., (IBerger et al.l l2007k ). For example, the offsets for GRBs 
060801, 061210, and 070429B are 19 ± 16 kpc, 11 + 10 kpc, 
and 40 + 48 kpc. I adopt 30 kpc as a typical upper limit on the 
offset for these 6 events. No follow-up observations are avail- 
able in the literature for most short GRBs with X-ray positions 
from 2008-2010. Finally, about 1/4 - 1/3 of all short GRBs 
discovered to date have only been detected in y-rays with po- 
sitional accuracies of a few arcminutes, thereby precluding a 
unique host galaxy association and an offset measurement. 

The cumulative distribution of p rojected physical o ffsets for 
the GRBs with HST observations dFong et al.l l2010l) . supple- 
mented by the bursts with offsets or limits based on optical af- 
terglow positions (070707, 070714B, 070724, 070809, 071227, 
080503, and 090510) is shown in Figure[T6l Also shown is the 
differential probability distribution, P(Sr)d(Sr), taking into ac- 
count the non-Gaussian errors on the radial offsets (see discus- 
sion in Appendix B of (IBloom et al.l 120021)) . The median for 
this sample is about 5 kpc (Fong et al., 2010). 

As evident from the preceding discussion, this is not a com- 
plete offset distribution; roughly an equal number of short 
GRBs have only limits or undetermined offsets due to their de- 
tection in just the X-rays or 7-ray£|. Taking these events into 
account, our most robust inferences about the offset distribution 
of short GRBs are as follows: 

• At least 25% of all short GRBs have projected physical 
offsets of < 10 kpc. 

• At least 5% of all short GRBs have projected physical off- 
sets of > 20 kpc. 

• At least 50% of all short GRBs have projected physical 
offsets of < 30 kpc; this value includes the upper limits 
for the hosts identified within XRT error circles. 

These robust constraints are marked in Figure [T6l 



^GRB 070714B is located at z = 0.923, while the redshift of GRB 070707 is 
not known. Based on the faintness of the host, R ^ 27.3 mag, we assume z = I 
to calculate the physical offset. 

^I do not consider the bursts that lack host searches since there is no a priori 
reason that these events (mainly from 2008-2010) should have a different offset 
distribution compared to the existing sample. 
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I next compare the short GRB offset distribution with th e off- 
sets of long GRBs from the sample of iBloom et alJ (120021) : see 
Figure [161 The offset distribution of long GRBs has been used 
to argue for a m assive star progenit or population, and against 
NS-NS binaries dBloom et al.ll2002h . The offset distribution of 
short GRBs is clearly shifted to larger physical scales. In par- 
ticular, the median offset for the long GRBs is 1.1 kpc, about a 
factor of 5 times smaller than the median value for short GRBs. 
Similarly, no long GRB offsets are larger than about 7 kpc, 
whereas at least some short GRBs appear to have offsets in ex- 
cess of 15 kpc. The significant difference between the offset 
distributions indicates that short GRBs do not arise from the 
same progenitor population as long GRBs. 

I further compare the observed distribution (and the ro- 
bust constraints outlined above) with predicted distributions 



for N S -NS binarie s in M il ky Way type galaxies (iBloom et al. 
19991: iFrver et al.1 1 19991: iBelczvnski etall |2006|), appropri 



ate for t he observed luminosities of short GRB host galaxies 
(IBergen 120091 There is good agree ment between the obse rved 
distribution and those predicted by iBloom et akl (Il999h and 



Belczynski et al.l (l2006h . The offset distribution of lFryer et al. 
(Il999h . with a median of about 7 kpc, predicts larger offsets 
and therefore provides a poorer fit to the observed distribution, 
which has a median of about 5 kpc. However, all three predicted 
distributions accommodate the offset constraints. In particular, 
they predict about 60 - 75% of the offsets to be < 10 kpc, about 
80 - 90% to be < 30 kpc, and about 10 - 25% of the offsets 
to be > 20 kpc. Thus, the projected physical offsets of short 
GRBs are consistent with population synthesis predictions for 
NS-NS binaries. However, the observations are also consis- 
tent with partial contribution from other progenitor systems for 
which kicks are not expected (magnetars, WD- WD binaries, 
accreting NS). 

5.3. Light Distribution Analysis 

In addit i on to the offset analysis in the previous section, 
Fong et akl (1201 Ol) studied the local environments of short GRBs 
using a comparison of the host brightness at the GRB location 
to the overall host light distribution. This approach is advan- 
tageous because it is independent of galaxy morphology, and 
does not suffe r from am biguity in the definition of the host cen- 
ter (see (Fru chter et al., 2006)). I note that for the regular mor- 
phologies of most short GRB hosts ( ^5.11) , the definition of the 
host center is generally robust, unlike in the case of long GRBs 
tFruch teretal., 2006; Wa inwright et a l., 2007). On the other 
hand, this approach has the downside that it requires precise 
pixel-scale positional accuracy. In the existing sample, this is 
only available for 6 short bursts. 

The fraction of total host light in pixels fainter than the after- 
gl ow pixel brightnes s for each host/filter combination is given 



Fong et akl (I2OIOI) . The cumulative light distribution his- 



in 

togram is shown in Figure [T71 The shaded histogram represents 
the range defined by the dual filters for 5 of the 6 bursts. The 
upper bound of the distribution is defined by the blue filters, 
indicating that short GRBs trace the rest- frame optical light of 
their hosts better than the rest-frame UV. This indicates that 
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Figure 17: Cumulative distribution of fractional flux at the location of short 
GRBs relative to their host light. For each burst the fraction of host light in 
pixels fainter than the GRB pixel location is measured. The shaded area is 
defined by the results for the two available filters for each short GRB. Also 
shown are data for long GRBs (dark gray line) and for cor e-collapse a n d Typ e 
la SNe (fight gray lines) from Fruchter et al. (2006) and JKellv et al.l i2008l) . 
The dashed line marks the expected distribution for objects which track their 
host light distribution. Short GRBs appear to under-represent their host light, 
while long GRBs tend to be concentrated in th e brightest regions of their hosts 
JFruchter et aD.l2006[) . From lFone et"S] ( l20ld) . 



short GRB progenitors are likely to be associated with a rela- 
tively old stellar population, rather than a young and UV bright 
population. 

The overall distribution has a median value of ^ 0.1 - 0.4 
(blue vs. red filters); namely, only in about one-quarter of the 
cases, 50% of the host light is in pixels fainter than at the GRB 
location. Thus, the overall distribution of short GRB locations 
under-represents the host galaxies light distribution, but traces 
the red light (old stars) more closely than the blue light (star 
formation). This is also true in comparison to the distribution 
for core-collapse SN e, which appear to track their host light 
(IFruchter et al.ll2006|). and even Type la SNe, which have a me- 
dian of about 0.4 (IKelly et al.Ll2008k . Thus, the progenitors of 
short GRBs appear to be more diffusely distributed than Type 
la SN progenitors. 



An extensive analysis of the brightness dist ribution at the 



location of long GRBs has been carried out by iFruchter et al 
(2006). These authors find that long GRBs are more concen- 
trated on the brightest regions of their hosts than expected from 
the light distribution of each host. In particular, they con- 
clude that the probability distribution of long GRB positions is 
roughly proportional to the surface brightness squared. As can 
be seen from Figure [T71 short GRBs have a significantly more 
diff'use distribution relative to the host light than long GRBs. In 
particular, for the latter, the median light fraction is about 0.85 
compared to about 0.25 ±0.15 for the short GRBs. 
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Figure 18: Images of the fields around the optical afterglow positions of short GRBs 061201, 070809, 080503, 090305A, and 090515 (top-left to bottom-right). 
These are the 5 bursts in the current sample that have optical afterglow positions and no coincident hosts to limits of > 25.5 mag. From .Berger. (2010.) . 



6. Is There Evidence for Large Progenitor Kicks? 

One of the predictions of the compact object coalescence 
model is that some binaries could merge at large separations 
from their host galaxies due to a combi nation of large kick 
velocities and a long merger timescale (JBloom et all 1 19991: 
Frver et all 1 19991: iBelczvnski eTall l2006h . With kicks of sev- 



eral hundred km s" and merger timescales of ~ Gyr, such a 
binary could travel several hundred kpc from its host, corre- 
sponding to ~ r at z ~ 1 and ~ 10' at z ~ 0.1. Such large 
offsets would not be expected in other progenitor scenarios. It 
is important to note, however, that if the typical kick velocities 
are < 10^ km s~\ an NS-NS/NS-BH system is likely to re- 
main bound to its host regardless of the merger timescale, and 
hence to reside at off sets of < tens of k pc (as already inferred 
for some short GRBs (iFong et alll2010h : FigureO. Similarly, 
short merger timescales (tens to hundreds of Myr) would also 
lead to relatively small offsets regardlss of the kick velocity. 

Clearly, an observational demonstration of a large offset is 
not trivial. Ideally, we would like to measure the redshift of 
the burst directly through afterglow spectroscopy and then as- 
sociate it with a galaxy at a large separation. However, to date, 
short GRB redshifts have been measured through their host as- 
sociations so this test is not possible. 

At a more tentative level, we can also investigate the large- 
scale environments around short GRBs that do not appear to 



spatially coincide with bright galaxy counterparts to assess the 
potential for a host with large offset. This is a particularly im- 
portant test if combined with the afterglow properties of short 
GRBs with and without coincident hosts. In the current sample, 
there are 5 cases of short GRBs with sub-arcsecond positions 
and no obvious bright hosto Below, I assess the possibility of 
large offsets for these bursts, and compare this with alternative 
explanations (e.g., a high redshift origin). The sample includes 
20 short GRBs with optical afterglows. Images of the fields 
around the 5 bursts with no coincident hosts are shown in Fig- 
ure [TSl hereafter, I denote these 5 bursts as Sample 2, with the 
remaining 15 bursts with coincident hosts designated as Sam- 
ple 1. The afterglow positions, as well as nearby galaxies with 
varying probabilities of chance coincidence are marked in Fig- 
ure [181 The limits at the positions of the afterglows range from 
> 25.5 to > 26.5 mag (Bergd,l2Q10|). 

In terms of the afterglow properties, the bursts in Sample 2 
have a median optical brightness that is 1 .4 mag fainter than the 
bursts in Sample 1 (Figure [191 and X-ray fluxes at 8 hours that 
are about a factor of 2 times smaller. The ir y-ray fluences are 
similarly smaller, by about a factor of 5 dBergeii l2010h . The 
differences in optical afterglow brightness can be due to lower 



^I do not discuss cases with only XRT positions of a few arcsec radius since 
those do not generally lead to significant off'sets and furthermore nearly always 
contain at least one possible host consistent with zero off'set. 
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Figure 19: Optical afterglow brightness on timescales of a few hours after the 
burst for short GRBs with detected afterglows (Sample 1: black squares; Sam- 
ple 2: red squares) or upper limits (gray triangles). The lines at the top right 
indicate the fading tracks for afterglow decay rates of of = -0.5 and -1. The 
right panel shows the projected histogram for the bursts with detected after- 
glows (hatched) and upper limits (open). The symbols mark the mean for each 
sample, and the vertical bar marks the standard deviation for Sample 1 . From 
TBerger i20ia) . 
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densities or higher redshifts for Sample 2. Both of these scenar- 
ios would also explain the lack of bright coincident hosts, since 
low densities may be indicative of large offsets and high red- 
shifts will lead to fainter host galaxies. The difference in X-ray 
brightness does not significantly constrain these two possibil- 
ities, while the fainter y-ray emission of the bursts in Sample 
2 points to high redshift as the likely explanation, since in the 
context of the standard GRB model the prompt emission flu- 
ence does not depend on the circumburst density. 

I next turn to an analysis of the large-scale environments of 
the bursts in Sample 2, particularly in comparison to the hosts 
of bursts in Sample 1 . None of the 5 bursts have coincident 
hosts to significantly deeper limits than the hosts in Sample 1 . 
I therefore investigate the possibility of large offsets through 
the calculation of chance coincidence probabilities for nearby 
galaxies, as well as the possibility of a high redshift origin. The 
chance coincidence probability for nearby galaxies depends on 
both their apparent magnitude and their distance from the opti- 
cal afterglow position. The expected number density of galax- 



01 galax- 
ies brighter than a measured magnitude, m , is (iHogg etaU I ????! : 
Bloom et al.l.l2002l:lBeckwith etal.l [20061) : 



cr(< m) = 



1 



0.33xln(10) 



>< lQ0.33(m-24)-2.44 ^^^^^^-2^ ^3) 



and the chance coincidence probability for a given separation, 
P(< SR), is then: 



P(< SR) = 1- e 



-7T(SRYo-(<m) 



(4) 



where for offsets substantially larger than the galaxy size, SR is 
the appropriate radius in Equation|4]( Bloom et al., 2002). 

The resulting distributions for each field are shown in Fig- 
ure [20l I include all galaxies that have probabilities of < 0.95. 



Figure 20: Probability of chance coincidence as a function of distance from 
a short GRB optical afterglow position for galaxies near the location of each 
burst. These are the galaxies marked in Figure [l^ In each panel I mark the 
galaxy with the lowest probability of chance detection with a circle. In 4 of 
the 5 cases, the lowest probability is associated with galaxies that are offset by 
.^^^lJ.^,1i^^^^^^^^' ^^^^ ^^^ nearest galaxies are offset by ^ 1.6-5.8''. From 
iBergeJ ( [20TqI) . 



I find that for 4 of the 5 bursts, faint galaxies (~ 25-26 
mag) can be identified within ^ 1.6 - 2'' of the afterglow 
positions, with associated chance coincidence probabilities of 
^ 0.1 - 0.2; in the case of GRB 090515 I do not detect any 
such faint galaxies within ^ 5'' of the afterglow position. For 
GRB 080503 I also include the galaxy at a n offset of 8^^ and 
mAB(F606W) = 27.3 ±0.2 mag identified by 'PerlexelalJ (|2009|) 
based on their deeper stack of HST observations. On the other 
hand, for 4 of the 5 bursts I find that the galaxies with the lowest 
probability of chance coincidence, ^ 0.03 - 0.15, are brighter 
objects with offsets of about 6 - 16'' from the burst positions; 
only in the case of GRB 080503 the lowes t chance coinciden ce 
is associated with the nearest galaxy (see (IPerley et al.[|2009k ). 
The use of a posteriori probabilities to assign unique galaxy 
associations is fraught with difficulties. First, for a given ap- 
parent brightness, galaxies located further away from the GRB 
position, potentially due to larger kicks and/or longer merger 
timescales in the NS-NS merger framework, have higher prob- 
abilities of chance coincidence. Since no a priori model- 
independent knowledge is available for the range of possible 
kicks and merger timescales, it is not possible to rule out galax- 
ies at very large offsets for which P{< SR) ~ 1. Indeed, a 



reasonable constraint of Vkick ^10 km s"^ and Tmerger ^10 
Gyr leads to only a weak constraint on the offset of < 10 Mpc. 
At z = O.l (z = 1) this corresponds to about 1.5° (0.3°), a pro- 
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jected distance at which nearly all galaxies will have a chance 
coincidence probability of order unity. However, this fact only 
means that we cannot rule out offsets that are even larger than 
inferred from the most likely host association. 

A second difficulty is that I am using angular offsets, which 
ignore the potential wide range of redshifts (and by extension 
also luminosities) of the various galaxies. For example, if the 
faint galaxies with small offsets are located at z ^ 1 , the cor- 
responding physical offsets are ~ 15 kpc, while if the galaxies 
at ~ 10'' offsets are located at z ~ 0.3, the offsets are only 
somewhat larger, ~ 30 kpc. A galaxy at an even lower red- 
shift, z ~ 0.1, with an offset of 50 kpc will be located about 
30'' from the GRB position and incur a large penalty in terms 
of chance coincidence probability. It is important to note, how- 
ever, that galaxies at lower redshift will generally have brighter 
apparent magnitudes, partially compensating for the larger an- 
gular separations (Equations [3] and |4l). In only a single case 
(GRB 070809) I find a galaxy with P(< SR) < OA Sit SR > 1' 
(which at z = 0.043 for this galaxy corresponds to a physical 
offset of about 100 kpc). 

A final complication, which is not unique to this subset of 
events, is that only projected offsets can be measured, SR = 
SR3J) X cos(0). The measured offsets can be used as lower limits 
on the actual offsets, while for the overall distribution it is pos- 
sible to apply an average correction factor of 7r/2, based on the 
expectation value for the projection factor, cos(^). 

Despite these caveats, it is possible to address the probability 
that all of the associations are spurious. This joint probability is 
simpl y the product of the individual probabilities (JBloomet al. . 
12002"). For the faint galaxies at small angular separations the 
probability that all are spurious is Paii ~ 8 x 10"^, while for 
the galaxies with the lowest probability of chance coincidence 
the joint probability is nearly 30 times lower, Paii ~ 3 x 10"^. 
Conversely, the probabilities that none of the associations are 
spurious are ^ 0.42 and ^ 0.59, respectively. These values indi- 
cate the some spurious coincidences may be present in Sample 
2. Indeed, the probabilities that 1, 2, or 3 associations are spu- 
rious are [0.40, 0.15, 0.027] and [0.34, 0.068, 0.006], respec- 
tively. These results indicate that for the faint galaxies it is not 
unlikely that 2-3 associations (out of 5) are spurious, while for 
the brighter galaxies 1-2 associations may be spurious. This 
analysis clearly demonstrates why a joint statistical study is su- 
perior to case-by-case attempts to associate short GRBs with 
galaxies at substantial offsets. 

Based on the possibility of association with the galaxies at 
separatio ns of - 10", I obtained redshift for three of these 
galaxies (JBergeii I2OIOI) . leading to a star forming galaxy 
at z = 0.111 (GRB 061201), an early-type galaxy at z = 
0.473 (GRB 070809), and an early-type galaxy at z = 0.403 
(GRB 090515). The fainter host at separations of a few arcsec 
likely reside at z ^ 1 . The redshifts provide an indication of the 
physical projected offsets ( ^6.21) . 

6. 1 . Undetected Faint Hosts at High Redshift ? 

The redshifts of the GRBs in Sample 2 can be constrained 
based on their detections in the optical band (i.e., the lack 
of strong suppression by the Lya forest). The afterglow of 



GRB 061201 was detected in the ultr aviolet by the 5w?WU V0T 
and it is therefore located at z ^ 1.7 (JRoming et al.Ll2006b . The 
remaining four bursts were detected in the optic al g- or r-band , 
and can therefore be placed at z ^ 3 or < 4.3 ; see lBergej (l2010h 
for details. 
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Figure 21: Main Panel: Host galaxy optical magnitude as a function of redshift 
for short GRB hosts (black squares), long GRB hosts (detections: gray circles; 
non-detections or no redshifts: gray triangles), and galaxies with a luminosity 
of 0. 1 - 1 L* (shaded region). The dashed lines mark the upper limits at the GRB 
positions for the short GRBs with no coincident hosts. The arrows mark the up- 
per limits on the redshifts of three bursts in Sample 1 with faint hosts, based on 
the detection of the afterglows in the optical band (i.e., lack of a strong Lyman 
break). If underlying host galaxies exist for Sample 2, their non-detection indi- 
cates z > 1.5 (for 0.1 L*) or > 3 (for L*). The alternative possibility that they 
are located at similar redshifts to the detected hosts, requires < 0.01 L*, but 
this does not naturally explain their fainter afterglows. Upper Panel: Projected 
redshift histogram for Sample 1 (black) and Sample 2 (dashed red limits) un- 
der the assumption that the hosts are 0.1 L* galaxies (z ~ 1.5) and L* galaxies 
(z ~ 3). From iBergen J201Q) . 

It is possible to place additional constraints on the redshifts of 
any underlying hosts using the existing sample short GRB host 
galaxies. In Figure [21] I plot the r-band magnitudes as a func- 
tion of redshift for all available short GRB hosts from Sample 
1. For the faint hosts without known redshifts (GRBs 060121, 
060313, and 070707), I place upper limits on the redshift us- 
ing optical detections of the afterglows. A wide range of host 
magnitudes, tab ~ 16.5-27.5 mag, is appare nt. I also plot th e 
r-band magnitudes of long GRB hosts ( Savaglio et al.L l2009b . 
as well as the r - z phase space that is traced by galaxies with 
luminosities of L = 0.1 - 1 L*. I find excellent correspondence 
between the hosts of long and short GRBs, and the phase-space 
traced by 0.1 - 1 L* galaxies, at least to z ~ 4. In the context 
of these distributions, the available limits for the short GRBs 
in Sample 2 translate to redshifts of z ^ 1.5 if they are 0.1 L* 
galaxies, or z ^ 3 if they are L* galaxies. The latter lower 
limits are comparable to the redshift upper limits inferred from 
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the afterglow detections. I note that for GRB Q8Q 5Q3, the Hni' 



its of ^28.5 mag and z ^ 3 from the afterglow I Perley et al 



I 1 - -- - ^ ' 

12009) place even more stringent limits on the luminosity of an 

underlying galaxy of < 0.1 L* galaxy. 

The possibility that the five bursts originated at z ^ 3 leads 
to a bimodal redshift distribution (Figure [21]). Nearly all of the 
bursts in Sample 1 with a known redshift (9/10) have z ^ 0.2- 1 , 
with a media n of (z) ^ 0.5; the sole exception is GRB 090426 
atz = 2.61 ( Ant onelH erall l2009l: ILevesque et all I2OIO). The 
three bursts with faint coincident hosts have upper limits of z ^ 
4 from afterglow detections, while lower limit of z ^ 1.5-2 can 
be placed on these hosts if they have L > 0.1 L*. Adding the 
Sample 2 bursts with the assumption that they have z ^ 3 will 
furthermore result in a population of short GRBs with a median 
of z ~ 3, and leave a substantial gap at z ~ 1-2 (Figure[2T]). If 
the 5 bursts are instead hosted by 0.1 L* galaxies, the inferred 
lower limits on the redshifts (z ^ 1.5) lead to a potentially more 
uniform redshift distribution. 

It is difficult to explain a bimodal redshift distribution with 
a single progenitor population such as NS-NS binaries, with- 
out appealing to, for example, a bimodal distribution of merger 
timescales. Another possibility is two distinct progenitor popu- 
lations, producing bursts of similar observed properties but with 
distinct redshift ranges. While these possibilities are difficult to 
exclude, they do not provide a natural explanation for the short 
GRB population. 

A final alternative explanation is that any underlying hosts 
reside at similar redshifts to the known hosts in Sample 1 
(z ~ 0.5), but have significantly lower luminosities of < 0.01 
L* . This scenario would not naturally explain why the bursts in 
Sample 2 have fainter optical and X-ray afterglows, as well as 
lower y-ray fluences. I therefore do not consider this possibility 
to be the likely explanation. 

6.2. Large Offsets? 

While higher redshifts may explain the lack of detected hosts, 
the fainter afterglows, and the smaller y-ray fluences of the 
bursts in Sample 2, this scenario sufl'ers from several difficulties 
outlined above. The alternative explanation is that the bursts oc- 
curred at significant off'sets relative to their hosts, and hence in 
lower density environments that would explain the faint after- 
glow emission (though possibly not the lower y-ray fluences). 
As demonstrated in the chance coincidence analysis, the ofl*- 
sets may be ~ T' (-15 kpc) if the bursts originated in the faint 
galaxies at the smallest angular separations, or ~ 10'' (~ 30-75 
kpc) if they originated in the brighter galaxies with the lowest 
probability of chance coincidence (Figure[20b. 

The projected physical ofl'sets for Sample 1 and Sample 2 
are shown in Figure [22l The mean and standard deviation for 
Sample 1 are SR = 4.2 ±3.8 kpc, and a log-normal fit results 
in a mean of log(SR) ^ 0.5 and a width of (Tiog(^/?) ~ 0.3. On 
the other hand, the bursts in Sample 2 have a mean off'set of 
about 19 kpc if they arise in the faint galaxies with small angular 
separation, or about 40 kpc if they arise in the brighter galaxies, 
pointing to distinct distributions. 

A similar result is obtained when considering the off'set nor- 
malized by each host's eff'ective radius. Re (as advocated by 
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Figure 22: Top: Histogram of projected physical offsets relative to the host 
galaxy center for short GRBs with coincident hosts (hatched black), and bursts 
with no coincident hosts if the galaxies with lowest chance coincidence proba- 
bility are the hosts (hatched red), or if the faint galaxies with smallest angular 
separation are hosts (open red); see Figure |20] The dashed line is a log-normal 
fit to the bursts with coincident hosts. Bottom: Same, but for offsets normalized 
relative to the host effective radii. Re. The dashed line is a log-normal fit to the 
bursts with coincident hosts. From lBergej ( |201Q() . 



Fong et al.l (I2OIOI) ). This quantity takes into account the vary- 
ing sizes of the hosts due to both intrinsic size variations and 
redshift eff'ects. It also gives a better indication of whether the 
burst coincides with the host light or is significantly off'set. As 
shown in Figure [22l the host-normalized off'sets of Sample 1 
have a mean and standard deviation of about 1 ± 0.6 7?^, and a 
range of about 0.2 - 2 7?^. A log-normal fit results in a mean of 
log(SR/Re) ~ and a width of cr\og(^sRiRe) ~ 0.2. The bursts 
in Sample 2 have much larger host-normalized off'sets, with 
(SR/Re) = 7.3 ± 2.3 if they originated in the galaxies with the 
lowest chance coincidence probability. Even if I associate the 
bursts with the nearest faint hosts, the distribution has a mean 
of about 4 Re, reflecting the fact that the efl'ective radii of the 
faint galaxies are smaller than those of the brighter ones. 
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Figure 23: Cumulative distributions of projected physical offsets for short 
GRBs with coincident hosts (black line), and combined with offsets for the 
hosts with the lowest probabilities of chance coincidence (thick red line) or the 
faint hosts with smallest angular offsets (thin red line). Also sho wn are pre- 
dicted distributions for NS-NS kicks from several models (Bloom e t all 1 199^ : 
iFrver et al , 1999; Belczvnski et al., 2006), and for dynamically-formed NS-NS 
binaries from globular clusters (shaded region marks a range of predictions for 
host galaxy masses of 5 x 10^^ - 10^^ M©; ( Salvaterra et al., 2010)). The models 
with kick velocities are in good agreement with the measured offset distribution 
for either set of galaxy associat ions, wh ile the globular clusters model provides 
a poor match to the data. From lBergen 12010.) . 



Thus, the distributions of physical, host-normaUzed, and an- 
gular offsets exhibit a clear bimodality when associating the 
bursts in Sample 2 with the galaxies at z ~ 0.1-0.5. The 
effect is still apparent, though less pronounced in the case of 
association with the faint galaxies at z ^ 1. Thus, if the offset 
scenario is correct, the resulting distributions point to a possi- 
ble bimodality rather than a single continuous distribution of 
offsets. 

The cumulative distributions of physical offsets for Sample 1 
alone, and in conjunction with the two possible offset groups for 
Sample 2 are shown in Figure |23l The combined distributions 
have a median of about 4 kpc, driven by the bursts with coinci- 
dent hosts. However, there is a clear extension to larger physical 
offsets in the case of association with the brighter galaxies, with 
about 20% of all objects having dR > 30 kpc. The cumulative 
distributions are particularly useful for comparison with NS-NS 
merger models since predictions exist for both the kick scenario 
and the globular cluster origin model. 



7. Implications for the Progenitors 

The extensive analysis of host galaxy properties and the sub- 
galactic environments of short GRBs presented above provides 
important insight into the nature of the progenitors. I address in 
particular the popular NS-NS merger model, as well as delayed 
magnetar formati on via WD-^Vp m ergers or WD accretion- 
induced collapse ( Metzger et al.ll2008i) . 



7.1. Host Galaxy Demographics and Properties 

The identified hosts of short GRBs with sub-arcsecond po- 
sitions (generally from optical afterglows) are dominated by 
star-forming galaxies with a ratio of about 4:1, although I note 
that the nature of the hosts of several short GRBs with sub- 
arcsecond positions remain unknown mainly due to their faint- 
ness and/or lack of deep imaging and spectroscopic observa- 
tions. Only if all the unidentified hosts are early-type galaxies, 
would we have a ratio of 1 : 1 . The putative hosts identified in co- 
incidence with XRT error circles exhibit a similar ratio of about 
4:1. This result is also supported by the m orphologica l analy - 
sis of short GRB hosts observed with HST dFongetaU l2010h . 
which indicates that the ratio of hosts with an exponential disk 
profile versus a de Vaucouleurs profile is 4:1. Thus, I conclude 
that the most conservative estimate of the ratio of star forming 
to elliptical hosts is about 1:1, but that if the well-studied (i.e., 
brighter) hosts are representative of the whole sample, than the 
ratio is about 4: 1 . A ratio of about 1 : 1 is expected if short GRBs 
select galaxies by stellar mass alone. Thus, the existing demo- 
graphics cannot rule out this scenario, but do appear to point to 
an over- abundance of short GRBs in late-type galaxies, possi- 
bly indicative of a short GRB rate that partially depends on star 
formation activity. 

Although star forming hosts appear to dominate the host 
population, it is clear from a comparison of their luminosities, 
metallicities, stellar masses, population ages, sizes, and star for- 
mation rates that they are distinct from the hosts of long GRBs. 
Namely, they exhibit lower star formation activity and appear 
to be larger, more massive, and dominated by more evolved 
stellar populations. This indicates that the short GRB progeni- 
tors are mostly related to the old stellar populations within their 
hosts, and perhaps only partially to the modest level of on-going 
star formation activity. This is further supported by the dearth 
of morphological galaxy merger signatures, which point to de- 
lays of > 1 Gyr relative to any merger-triggered star formation 
episodes. These conclusions are in direct contrast to the mas- 
sive star progenitors of long GRBs, which select galaxies by 
star formation (and perhaps metallicity). While the relation to 
old stellar populations does not rule out models such as WD- 
WD mergers or WD/NS AIC, it does disfavor young magnetars. 

As an aside, I note that this result demonstrates that caution 
should be taken with the proposed re-classification of short and 
long GRBs into Type I an d II events, marking old a nd young 
progenitors, respectively rZhan g et al.L 120071 l2009k . Such a 
new bimodal classification may lead to the erroneous conclu- 
sion that short GRBs in late-type galaxies (even if they track 
the on-going star format ion activity) share the same progeni- 
tors as long GRBs (e.g., (Ivirgili et al.[ l2009b ) since both would 
be classified as Type II. At the very least, such a new classi- 
fication scheme may require a further breakdown of the Type 
II events into those that result from massive stars versus those 
that simply track star formation activity with a modest delay, 
e.g.. Type Ila and lib. Clearly, this is beyond the scope of the 
current short GRB sample. 

In the context of the NS-NS/NS-BH merger models, the 
host demographics, coupled with the redshift distribution. 
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prov ide a rough con s traint on the typical me r ger time de 



lay dGuetta and PiranL l2QQ6l: iNakar et all [2OO6I: Berger et al 



2007; 



Zheng and R amirez-Ruiz!, l2007h . In the formulation of 



Zheng and Ramirez-Ruiz (2007), with a merger timescale prob- 
ability distribution of P(t) oc t", the dominance of star forming 
hosts with a ratio of about 4:1 indicates ^ < -1. I note that 
this is only a rough estimate since the calculation is appropri- 
ate for z ~ 0, while the host redshift distribution extends at 
least to z ~ 1 . Similarly, a comparison of the observed red- 
shift distribution (with a range of z - 0.2 - 2 and a median 
of ~ 0.7) to the models of iNakar et all (l2006h (their Figure 2), 
indicates a typical progenitor merger timescale of about 3 Gyr 
(for a log-normal m erger timescale di s tributi on). A compari- 
son to the models of lOuetta and PiranI (l2006h (their Figure 4), 

indicates a value of n 1 (for a power law merger timescale 

distribution), in good agreement with the results from host de- 
mographics. 

Thus, I conclude that based on the nature of short GRB host 
galaxies and their redshift distribution, the merger timescale 
distribution for the NS-NS/NS-BH coalescence models can be 

characterized by a power law index of 1 (power law model) 

or a typical value of ~ 3 Gyr (log-normal model). If the subset 
of short GRBs in late-type galaxies tend to track star forma- 
tion activity, their typical p rogenitor ages are instead -0.2 Gyr 
( Leibler and BergeA l2010h . It should be noted that the merger 
timescale distribution for the admittedly small sam ple of Galac- 
tic NS-NS binar ies has been claimed to be 1 (IPiranL 1 19921 : 

Champion eTaP. 12004) . 



7.2. Short GRB Offsets 

The differential offsets measured from the HST observations 
provide the most precise values to date for short GRBs, with a 
total uncertainty of only ~ 10-60 mas (< 1 pixel), correspond- 
ing to ~ 30 - 500 pc. None of the offsets are smaller than ~ 1 
kpc, while this is the median offset for long GRBs. On the other 
hand, a substantial fraction of the offsets are only a few kpc, in- 
dicating that at least some short GRBs explode within the stel- 
lar component of their hosts (rather than their extended halos). 
The median offset from the HST observations supplemented by 
ground-based data is about 5 kpc (Figure [161), roughly 5 times 
larger than for long GRBs. 

As discussed above, the observed offset distribution is in- 
complete since about 1/4 - 1/3 of all short GRBs have only 
y-ray positions (~ 1 - 3'), and a similar fraction have only XRT 
positions, which generally lead to a range of ~ 0-30 kpc. 
Taking these limitations into account I find that the most robust 
constraints on the offset distribution are that > 25% of all short 
GRBs have offsets of < 10 kpc, and that ^5% have offsets 
of > 20 kpc. In addition, for the current sample of of short 
GRBs with sub-arcsecond afterglow positions and no coinci- 
dent bright hosts, I find evidence for offsets of ~ 50 kpc (® . 

Both the observed offset distribution and the various con- 
straints are in good agreement with predictions for the off- 
set distribution of NS-NS b inaries in Mil ky Way type galaxie s 
JBloom et al.lll999l:lFrver~et al, 1999: Belczvnski et al.ll2006h . 



young magnetar flares, cannot be ruled out. The existence of 
modest (-10 kpc) and perhaps large (~ 50 kpc) offsets in the 
sample suggests that these latter models are not likely to ac- 
count for all short GRBs. 

In the context of implications for the pro genitor population , 
the study of short GRB physical offsets by iTroja et al.l (l2008h 
led to the claim that short GRBs with extended X-ray emis- 
sion have systematically smaller offsets than those with only 
a prompt spike, possibly due to a systematic difference in the 
progenitors. The HST sample of iFong et al.l (.2010.) includes 
three short GRBs with significant extended emission (050709, 
050724, and 061006), and one burst (060121) with possib le 
extended emission (4.5cr significance; (JDonaghy et al.Ll2006b ). 
The physical offsets of these bursts are about 3.7, 2.7, 1.3, and 1 
kpc, respectively, leading to a mean offset of about 2.2 kpc. The 
physical offsets of the bursts without extended emission, but 
with precise afterglow positions (051221, 060313, and 061201) 
are 2.0, 2.3, and 14.2 or 32.5 kpc, respectively. The two events 
with no extended emission and with XRT positions (050509B 
and 051210) have offsets of about 54+12 and 28 ± 23 kpc, 
respectively. Including the ground-based sample with optical 
afterglow positions, the bursts with apparent e xtended emis 
sion ( 070714B, 071227, 080513, and 0905 10; dBarbier et al 



2007l:ISakamoto et al.l.[2QQ7l:IUkwatta et al. 1120091: IPerlev et al. 



However, at the present a partial contribution from other pro- 
genitor systems, such as delayed magnetar formation and even 



20091) ) have offsets of < 4, 14.8, ~ 20, and ~ 5.5 kpc, while the 
bursts without extended emission (070724 and 070809) have 
offsets of 4.8 and ~ 6.5 kpc. Thus, in the sample of events 
with sub-arcsecond positions, 6/8 bursts with extended emis- 
sion have offsets of < 5 kpc and 2/8 have likely offsets of 
~ 15-20 kpc. In the sample without extended emission, 4/5 
have offsets of < 6 kpc and 1/5 has a likely offset of ~ 14-32 
kpc. Thus, at the present it does not appear that there is a sig- 
nificant difference in the two offset distributions. 

The inclusion of events with only XRT positions does not al- 
ter this conclusion. In particular, of the subset with no extended 
emission only GRB 050509B is likely to have a significant off- 
set, while GRBs 051210, 060801, and 070429B have off'sets 
(28 ± 23, 19 ± 16, and 40 ± 48 kpc, respectively) that are con- 
sistent with zero. Similarly, GRB 061210 with extended emis- 
sion has an offset of 1 1 ± 10 kpc. A continued investigation of 
the difference between short GRBs with and without extended 
emission will greatly benefit from the use of host-normalized 
offsets, which take into account the individual hosts' effective 
radii. 

7.3. Large Offsets? 

The sample of short GRBs with optical afterglows represents 
about 1/3 of all short bursts, and may thus not be fully repre- 
sentative. One often-discussed bias is that the bursts with op- 
tical afterglows require a high circumburst density, and there- 
fore have negligible offsets. However, from the analysis in ^ 
(Berger, 2010) it is clear that one explanation for the lack of co- 
incident hosts for the bursts in Sample 2 is indeed large offsets, 
despite their detection in the optical band. 

As shown in Figure[23l the NS-NS merger model predictions 
have a median of about 6 kpc, compared to about 4 kpc for 
the observed sample. On the other hand, the models predict 
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10 - 20% of offsets to be > 30 kpc, in good agreement with 
the observed distribution in both the ~ 15 kpc and ~ 40 kpc 
scenarios (©.I note that the overall smaller offsets measured 
from the data may be due to projections effects. Indeed, the 
mean correction factor of njl nicely reconciles the theoretical 
and observed distributions. 

In ^ I noted a bimodality in the physical and host- 
normalized offsets for Sample 1 and Sample 2 (Figure |22]). In 
the framework of NS-NS binary kicks this bimodality may indi- 
cate that the binaries generally remain bound to their host galax- 
ies, thereby spending most of their time at the maximal distance 
defined by Jmax = 2GMhost/v?;^ T. (i-^-. with their kinetic energy 
stored as potential energy; (.Bloom et al.l l2007h ). This would 
require typical kick velocities of less than a few hundred km 
s-i. 

I further compare the observed offset distribution to predic- 
tions for dynamically formed NS-NS binaries in globular clus- 
ters, with a range of host galaxy virial masses of 5 x 10^^ - 10^^ 
M© (.Salvaterra et al... .2010.) . These models predict a range of 
only ^ 5 - 40% of all NS-NS mergers to occur within 10 kpc 
of the host center, in contrast to the observed distribution with 
about 70% with dR<\Q kpc. I stress that this result is indepen- 
dent of what offsets are assigned to the bursts in Sample 2 since 
they account for only 1/4 of the bursts with optical afterglows. 
On the other hand, the globular cluster origin may account for 
the bimodality in the physical and host-normalized offsets (Fig- 
ure [22|, with the objects in Sample 2 arising in globular clusters 
and the objects with coincident hosts arising from primordial 
NS-NS binaries. This possibility also agrees with the predicted 
fraction of dynamicall y-formed NS-NS binaries of- 10-30% 
( Grindlay et al.l l2006b . The cumulative offset distributions for 
Sample 2 alone (assuming the hosts are the galaxies with the 
lowest probability of chance coincidence) is well-matched by 
the range of predictions for dynamically-formed NS-NS bina- 
ries in globular clusters (Figure [23]). In this scenario, however, 
the implication is that short GRBs outside of globular clusters 
do not experience kicks as expected for NS-NS binaries since 
the largest measured offset is only 15 kpc. 

Unless the populations of short GRBs with only X-ray or y- 
ray positions have fundamentally different offset distributions, 
I conclude that the measured offsets of short GRBs and the pre- 
dicted offsets for NS-NS kicks are in good agreement, if when 
treating all short GRBs with optical afterglows as a single pop- 
ulation. Alternatively, it is possible that the bimodal distribu- 
tions of physical and host-normalized offsets point to a progen- 
itor bimodality, with the bursts in Sample 2 originating in glob- 
ular clusters. 



7.4. Relation to the Host Galaxy Light Distribution 

In addition to the projected offsets relative to their host cen- 
ters, it is apparent that short GRBs are more diffusely dis- 
tributed relative to their host light than long GRBs. In partic- 
ular, the locations of short GRBs under-represent their over- 
all host light distributions, even in comparison to core-collapse 
and Type la SNe. On the other hand, it appears from the cur- 
rent small sample that short GRBs are better tracers of their 



hosts' rest- frame optical light than UV light. This result indi- 
cates that short GRBs arise in locations within their hosts that 
trace the distribution of older stellar populations, and clearly do 
not trace the sites of active star formation. This result provides 
strong support to the claim that although most short GRB hosts 
are star forming galaxies, the bursts themselves are not related 
to the star formation activity (©. 

At the present, the sample of events with sufficiently precise 
astrometry to determine the burst locations at the level of < 1 
HST pixel is very small (6 events). It is therefore not possi- 
ble to draw conclusions about the fraction of short GRBs that 
are associated with old stellar populations as opposed to young 
populations (as expected for young magnetars). Luckily, there 
are at least 10 additional events for which these measurements 
can be made with future HST observations. At the present, I 
conclude that the stronger correlation of short GRBs with the 
rest- frame optical light than UV light of their hosts is indicative 
of a dominant old progenitor population. 



8. Conclusions 

While the sample of short GRBs with afterglow positions 
is still significantly smaller than that of long GRBs, we have 
made significant progress in understanding their galactic and 
sub-galactic environments. The results of host galaxy imag- 
ing and spectroscopy, including high-resolution imaging with 
HST, point to an association of short GRBs with old stellar 
populations within a range of normal star forming and elliptical 
galaxies. In nearly every respect (star formation rates, metal- 
licities, sizes, offsets, light distribution) the environments of 
short GRBs are distinct from those of long GRBs, indicating 
that their progenitors are not related to a young progenitor pop- 
ulation. 

As I showed through the study of short GRB offsets, host 
galaxy demographics, and the redshift distribution, the current 
observations are fully consistent with NS-NS/NS-BH binary 
mergers. However, a partial contribution from other (mainly 
old) progenitor channels (e.g., WD- WD mergers leading to 
magnetar formation, WD/NS AIC) cannot be ruled out at the 
present. Currently, we do not have conclusive evidence for sig- 
nificant progenitor kicks, which are only expected in the coa- 
lescence model. Still, a few events with sub-arcsecond optical 
positions do not directly coincide with bright host galaxies, and 
yet reside within tens of kpc from bright, low-redshift galaxies. 
This may be suggestive of progenitor kicks, but it is also possi- 
ble that these bursts are associated with fainter hosts (likely at 
higher redshift) with marginal offsets. 

With continued vigilance, and a short GRB discovery rate 
of about 1 event per month, we are likely to gain further in- 
sight into the nature of short GRB progenitors in the next few 
years, possibly with the first detections (or significant limits) of 
gravitational waves. As argued in this review, host galaxy ob- 
servations of existing and future events will play a central role 
in our on-going quest to determine the identity of short GRB 
progenitors. 
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